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Abstract 

We investigate the replica trick for the microscopic spectral density, Ps{x), of the 
Euclidean QCD Dirac operator. Our starting point is the low-energy limit of the QCD 
partition function for n fermionic flavors (or replicas) in the sector of topological charge 
u. In the domain of the smallest eigenvalues, this partition function is simply given by a 
U{n) unitary matrix integral. We show that the asymptotic behavior of Ps(x) for x ^ oo 
is obtained from the n ^ limit of this integral. The smooth contributions to this series 
are obtained from an expansion about the replica symmetric saddle-point, whereas the 
oscillatory terms follow from an expansion about a saddle-point that breaks the replica 
symmetry. For z/ = we recover the small-x logarithmic singularity of the resolvent by 
means of the replica trick. For half integer z/, when the saddle point expansion of the U{n) 
integral terminates, the replica trick reproduces the exact analytical result. In all other 
cases only an asymptotic series that does not uniquely determine the microscopic spectral 
density is obtained. We argue that bosonic replicas fail to reproduce the microscopic 
spectral density. In all cases, the exact answer is obtained naturally by means of the 
supersymmetric method. 
PACS: ll.30.Rd, 12.39.Fe, 12.38.Lg, 71.30.+h 

Keywords :Iieplica Trick; QCD Dirac Spectrum; Chiral Random Matrix Theory.. 



^On leave from UNESP (Guaratingueta-Brazil) 



1 Introduction 



In disordered systems, the ensemble average of the logarithm of the partition function 
cannot be evaluated directly in most cases. Two widely used methods, the replica trick 
^ and the supersymmetric method 0, ^ §], have been proposed to circumvent this 
problem. In the replica trick the ensemble average of the logarithm of the partition 
function is written as 

Vn _ 1 

i.e., logZ is calculated by the analytic continuation of the n-dependence of n replicated 
partition functions to n 0. If Z is given by a fermion determinant, can be written 
as an integral over n replicated Grassmann fields, and therefore, this method is known as 
the fermionic replica trick. Alternatively, log Z can be written as 

1 — 7^" 

(logZ) = lim( ). (2) 

In this determintal Z can be expressed as an integral over n replicated complex 

fields and this limit referred to as the bosonic replica trick. Both for the fermionic and 
for the bosonic replica trick, the average partition function can mapped onto a non-linear 
cr-model which is amenable to a saddle-point expansion. 

In the supersymmetric method the disorder average is performed for the ratio 



{dj\ogZ{J))j^,= d. 



I Z{J) ' 
\Z(J = 0), 



(3) 

j=o 



If Z{J) is given by a determinant, the numerator can be expressed as a fermionic in- 
tegral whereas the denominator can be written as a bosonic integral. For J = the 
partition function is thus invariant with respect to superunitary transformations that mix 
the fermionic and bosonic fields. Based on this symmetry, the supersymmetric partition 
function can be mapped onto a supersymmetric non-linear a— model p, Q which can be 
used to derive exact analytical expressions for spectral correlation functions. 

he replica trick has two obvious advantages. First, it can be applied to cases where Z 
cannot be expressed as a determinant (as for example is the case in the theory of spin- 
glasses [^), and, second, it is possible to calculate the logarithm of the partition function. 
One disadvantage of this method is that, in order to take the replica limit, the n depen- 
dence of the replicated partition function has to be known analytically. Therefore, the 
application of replica trick is limited to perturbative expansions of the partition function. 
A much more serious problem of this method is that the continuation of the n-dependence 
to n — i> is not unique. For example, a term of the form ~ sin nvr contributing to Z^ van- 
ishes for integer n but gives rise to a nonzero result in the replica limit. The failure of the 
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replica trick was first noticed in the theory of spin glasses where a replica symmetric 
minimum of the free energy resulted in a negative entropy 0. However, in that case the 
problems could be resolved by means of an elaborate scheme of replica symmetry break- 
ing 0. More recently, the replica trick was criticized because of its failure to reproduce 
the oscillatory contributions to random matrix correlation functions 0. Notice however 
that the non-oscillatory contributions to the two-point function were reproduced correctly 
P, P . Another example for which the replica trick may be problematic are nonhermitian 
Random Matrix Theories with eigenvalues scattered in the complex plane |T^, but we 
will not study such theories in this article. 

The advantage of the supersymmetric method is that it is possible to derive non- 
perturbative analytical results. This has been shown convincingly for the calculation of 
spectral correlation functions of random matrix-ensembles ^, |^, |ll], |12| . A disadvantage 
is that it requires some familiarity with supermathematics, but for perturbative expansions 
this method is no more complicated than the replica trick. In the example that will be 
discussed in this article, the exact super-symmetric calculation is actually much simpler 
than the perturbative replica calculation. A second disadvantage is, that because the 
average partition function is normalized to unity, one does not have access to the average 
free energy. 



Recently, the replica trick was revived in an article by Kamenev and Mezard |13 



They found that, in order to reproduce the oscillatory terms in spectral two-point corre- 
lation function, saddle points with broken replica symmetry had to be taken into account 
T3[| . For the Gaussian Unitary Ensemble they found the exact analytical result. As ex- 



plained in an article by Zirnbauer , the reason for this miracle is a consequence of the 
Duistermaat-Heckman theorem |T^ which is applicable to the a— model for the two-point 
function of the Gaussian Unitary Ensemble. This theorem states the conditions under 
which an integral is localized on its critical points so that a saddle-point approximation 
becomes exact. For the Gaussian Orthogonal Ensemble and the Gaussian Symplectic 
Ensemble the replica trick could only reproduce the asymptotic expansion of the two- 
point spectral correlation function for large energy differences [|T6i [T^ . We remind the 
reader that, unless we know the analytical properties of a function in the complex plane, 
it cannot be reproduced from its asymptotic series. A version of the replica trick that does 
not rely on the non-linear a-model but instead on orthogonal polynomials was shown to 



reproduce the exact correlation functions for all Gaussian ensembles |]18[. In this article 
we will not discuss this variant of the replica trick which is not an alternative to the 
orthogonal polynomial method. 

To investigate the replica trick we have chosen the microscopic spectral density of 
the QCD Dirac operator which is defined as the spectral density near zero on the scale 
of the average level spacing. The reason is three-fold, i) The ra-fold replicated partition 



3 



function is the QCD partition function with n flavors. Its low-energy limit, relevant for the 
microscopic spectral density, is known analytically. Because of spontaneous breaking of 
chiral symmetry in QCD, it is given by a partition function of weakly interacting Goldstone 
bosons and, on the scale of the average level spacing, it can be reduced to a unitary matrix 
integral which can be evaluated analytically for any number of flavors [|19|, ^ This, 
so called finite volume chiral partition function has been investigated in great detail, also 



in the context of one-link integrals in lattice QCD [|T9|, |2^. ii) Because the eigenvalues 
occur in pairs ±A, the level repulsion of the eigenvalues leads to a nontrivial oscillatory 
behavior of the microscopic spectral density, iii) The low energy partition function can be 
derived as a function of two integer valued parameters, the topological charge u and the 
number of physical flavors Nf (in this article only Nf = will be discussed) and can be 
trivially continued to non-integer u. For half-integer i>, the saddle point expansion of the 
U (n)-integrals terminates for finite positive integer values of n. This is closely related to 
Duistermaat-Heckman localization, where the leading order saddle point approximation 
is exact such as for the cr— model with n fermionic replicas of the two-point function of 
the Gaussian Unitary Ensemble [l^. Based on recent work |14[, we expect that the 
replica trick with replica-symmetry breaking gives the exact result in this case. 

The replica limit of the finite volume chiral partition function was first studied in 



24^ . It was found that the asymptotic expansion of the valence quark mass dependence 



of the chiral condensate in the microscopic region (i.e. for large valence quark masses in 
units of the average eigenvalue spacing) was reproduced by the replica trick. The small 
mass expansion was obtained up to logarithmic singularities which are essential for the 
calculation of the spectral density. Recently, the partially quenched supersymmetric chiral 
Lagrangian EB, [Tl|, IT2|, was formulated in terms of the replica trick [p5|]. 



In this article we analyze the oscillatory contributions to the microscopic spectral 
density in the framework of the replica trick. In section 2 we discuss the chiral sym- 
metries for bosonic replicas and introduce the supersymmetric quenched low energy chi- 
ral partition function. The calculation of the microscopic spectral density by means of 
the supersymmetric method is given in section 3. This calculation illustrates that the 
compact/non-compact structure of the final result for the resolvent appears naturally in 
the supersymmetric method. In section 4.1 we derive the large-mass asymptotic behavior 
of the microscopic spectral density by means of the fermionic replica trick. The small- 
mass behavior is discussed in section 4.2. Bosonic replicas are discussed in section 5 and 
concluding remarks are made in section 6. In the Appendix we the derive trace correlators 
necessary for the replica calculation to fourth order in the inverse mass. 
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2 Chiral Symmetry 



In this section we discuss chiral symmetry for bosonic and fermionic rephcas, as well as for 
the super symmetric partition function. To illustrate the difference in flavor symmetries 
between bosonic and fermionic replicas we present a detailed discussion for the case of 
one flavor or replica. 

The question we wish to address in this article is whether the spectrum of the QCD 
Dirac operator in the sector of topological charge v can be obtained from the QCD 
partition function with n additional replica flavors with quark mass z. This partition 
function is given by 

Zl^f+^){z) = f[dA], deC{tp + z)Yl det{tl/) + mj) e^^^*^'^] , (4) 
J f=i 

where mi, ■ ■ ■ ,171^^. are the usual quark masses. The integral is over all gauge fields in 
the sector of topological charge u (which is chosen positive in this article) and is weighted 
by the Yang-Mills action. The replica limit of the resolvent or the chiral condensate is 
defined by 

= -^Tr / = lim ll|-lnZf . (5) 

Here, V4 is the Euclidean 4- volume. The spectral density, which in terms of the eigenvalues 
i\k of the Dirac operator is defined by 

p{X)=Y.6{X-Xk), (6) 

k 

follows from the discontinuity across the imaginary axis 

^ = i-[S(zA + e)-S(zA-e)]. (7) 
Below we only consider the dimensionless ratio 

Ew 1 r,,.p(A) 



V4S0 ^-00 z + iX 
°° 1 u 1 

-00 V42^Q I/4L0 V^I^qZ + lU 

as a function of the microscopic variable I^Sq-z. The chiral condensate, Eq, is defined as 
the limit of for z close to z = but many level spacings away from the center of 
the spectrum. As one can see from the second equality, it can be expressed as an integral 
over the microscopic spectral density defined by 

Ps{u) = -^p{-^). (9) 
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For spontaneously broken chiral symmetry, the spacing of the eigenvalues is given by 
7r/EoV4 so that @ is stable in V4 and can be calculated in the thermodynamic limit. 

Because the resolvent has a cut along the imaginary axis it is sometimes more 
convenient to use the relation G{iX — e) = —G{~iX + e) to rewrite the discontinuity as 

^ = i^[S(a + 6) + SHA + 6)]. (10) 

We thus only need to calculate the resolvent in the half-plane Re(z) > 0. If the argument 
represents the microscopic variable, this relation gives us the microscopic spectral 
density. 

The reason of working with the partition function (^) is that in the phase of sponta- 
neously broken symmetry its low energy limit in entirely determined by chiral symmetry 
and is a partition function of weakly interacting Goldstone modes (or pions). We are 
interested in the kinematical domain ||20| 

> Vl^^ > 1/Aqcd. (11) 

Because V^^^ 1/Aqcd, only the Goldstone modes contribute to the mass- dependence 
of the partition function . For quark masses for which the Compton wavelength of the 
Goldstone modes is much larger than the size of the box ( l/m^r ^ V^^^), the kinetic term 
of the chiral Lagrangian can be ignored, and only the constant fields contribute to the 
mass dependence of the low-energy partition function. Therefore, as we will see next, in 
the domain ([TT|) the QCD parition function can be reduced to a unitary matrix integral. 

In QCD with n fermionic flavors the chiral symmetry group is given by Uv{n) x UA{n). 
A Ua{^) subgroup of the axial symmetry group is broken by the anomaly. The remain- 
ing axial symmetry group is broken spontaneously by the formation of a nonzero chiral 
condensate and the vector symmetry group, Uy{n), remains unbroken. The Goldstone 
manifold is thus given by the axial group SUA{n). The low energy limit of the QCD par- 
tition function is uniquely fixed from the requirement that its transformation properties 
under the chiral symmetry group are the same as for full QCD. Taking into account the 
anomaly, one finds, for n flavors all with mass m, in the sector of topological charge u, 



the low-energy flnite volume partition function [20 



Zl^\x) = I dU (detf/)^ei^(^+^"), (12) 

JU£U{n) 

where x = ml^So. This partition function is valid in the range ([Tl|). 

In the domain ([TT|) , the low energy partition function can also be obtained from 
a Random Matrix Theory with the global symmetries of the QCD partition function 
. In this theory, the matrix elements of the Euclidean Dirac operator are replaced by 
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independently distributed Gaussian random variables. In the sector of topological charge 
z/, the Dirac matrix is thus has given by 



D = . " , " . (13) 




where W is a.n n x [n + u) matrix, and the integration over the gauge fields is replaced by 
an integration over the probability distribution of the matrix elements. If W is complex 
and the probability distribution is a function of traces Tr(Vr^VF)^, this ensemble is known 
as the chiral Unitary Ensemble (chUE) or the chiral Gaussian Unitary Ensemble (chGUE) 
if the distribution of the matrix elements is Gaussian. As was shown in , the statistical 



properties of the smallest eigenvalues of D do not depend on the details of the probability 
distribution of the matrix elements. 

In order to analyze the chiral symmetry for bosonic replicas let us first discuss the flavor 
symmetries for one flavor. The fermion determinant that occurs in the QCD partition 
function can be written as a integral over Grassmann variables. 

det{D + m)= J dxdxe- / '^'^x(i5+n^)x _ (14) 



- / d*x<P*{D+m)(l> (ig) 



The inverse determinant can be written as an integral over bosonic integration variables 

1 

det(L' + m) " 

In a well-defined theory the functional integral has to be convergent. This is automat- 
ically the case for the Grassmann integration in (n), but the bosonic integrals in ([T5|) 
is only convergent for positive m. The symmetries of the partition function should be 
compatible with these convergence requirements. In particular, 0* should be identified 
with the complex conjugate of (p, and not as an independent integration variable such as 
the fermionic variables x ^^id x- 

If we decompose the spinors according to the block structure of the Dirac operator 
the vector symmetry of the massive theory with one fermionic replica (0) is given by 





(16) 

and the axial U{1) symmetry of the massless theory can be written as 

(17) 





The U{1) vector symmetry of the massive bosonic theory (|T5D is the same as for the 
fermionic theory 

(t]^e-(tA, ftiUe--ftiV (18) 
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This transformation does not affect the complex conjugation properties of 0. However, the 
axial transformation (17) applied to the bosonic fields affects their complex conjugation 
properties. In this case the axial transformation that is compatible with the convergence 
of the bosonic integral is given by 

(19) 

The axial symmetry group is therefore not f/(l) but instead GI{1) /U{1). Of course, this 
axial transformation is also a symmetry of the fermionic partition function. 

For n bosonic flavors the vector symmetry is U{n) whereas the axial symmetry is 
given by 

(20) 






^1 VI 

^2 y I 



with the matrix H containing only real elements. The axial symmetry group is thus given 
by the coset Gl{n)/U{n). An explicit parameterization of this coset is given by AA'^ with 
A G Gl{n). We expect that the axial symmetry of the bosonic partition function 

is broken in the same way as in the fermionic case with a G/(l)/?7(l) coset broken explicitly 
by the anomaly and the remaining part of the coset broken spontaneously by the chiral 
condensate. In absence of explicit symmetry breaking, the Goldstone manifold is thus 
given by Gl{n)/U{n). 

The mass term introduced according to 

n 

Y: 0t'Mfc,0l + 0fM,V^ (22) 

k,l=l 

is invariant under the axial transformation ( pOD provided that the mass matrix is trans- 
formed at the same time as 

M e-^Me-^, e^M^e". (23) 

The bosonic partition function in the sector of topological charge v transforms covariantly 
with dH) 

^ det(e2^^)Z(-"\ (24) 

The low energy limit of bosonic partition function is uniquely fixed by the requirement 
that its transformation properties are the same as of Z^^^K In the sector of topological 
charge u it is given by 

= / dU (det Uy e^Tr(MU+MW-^) (25) 

JU£Glin)/U{n) 
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The measure dU is the invariant Haar measure. Below we consider only the case with a 
diagonal mass matrix with all nonzero matrix elements equal to m and use the definition 
X = mV^'Eo- This partition function is valid in the kinematical domain (0) where constant 
Goldstone fields are the only relevant degrees of freedom. 

Since fermionic integrals are always convergent the partition function is invariant under 
both compact {U{n)) and non-compact {Gl{n) /U{n)) axial transformations. However, 
the small mass behavior of the non-compact partition function is singular because of 
the volume of the non-compact group diverges. This is not the case for the fermionic 
partition function and therefore this Gl{n)/U{n) is not an admissible parameterization of 
the Goldstone manifold. The correct parameterization is given by the compact manifold 
U{n). At a more technical level, this follows from the fact that the transformations that 
lead to the non-compact integral are only legitimate for bosonic quarks. For fermionic 
quarks one necessarily finds a compact effective partition function. 

In the supersymmetric method the generating function of the quenched resolvent is 
given by 

/ det(ip + z + J) \ 



It can be written as a superintegral 
Z^{M, M^) = { f d(f)d(f)*dxdx exp 




Both M and are given by diag(z + J,z), but in order to study the transformation 
properties of the partition function, we keep them as general matrices. 

With spontaneously broken axial symmetry, the bosonic part of the Goldstone mani- 
fold is U{1) X G/(l)/f/(l). Because, the partition function is invariant under super-unitary 
transformations, the full symmetry group is given by the maximum Riemannian subman- 



ifold of GZ(1|1) ig. In], |T1 (we will denote this manifold by This manifold can 

be parameterized as 

U=('^l Z)- (28) 



For zero topological charge the generating function (p6D is invariant under Glii{l\l) x 
if at the same time the mass matrix is transformed as 

M U^^MUl, ^ Ul^M^Ur. (29) 

For nonzero values of v the generating function (^6[) is not invariant under (|29|) but 
transforms according to 

Z,{z, J) Sdet"(f/^if/i)Z,(2, J). (30) 



(27) 
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The low-energy partition function is obtained from the requirement that it should have 
the same transformation properties as the QCD partition function (pGl). In the sector of 
topological charge i/, it is given by 

Z,{z,J)= f df/Sdet'^(f/)e^str(MC/+Mtt/-)_ ^g^^ 

JueGi{i\i) 

The integration is over the Haar measure of Below, we only the consider the case 

of a diagonal mass matrix with both M and M"^ equal to diag(z + J, z). 

An amusing observation is that the topological charge in QCD partition function is 
discrete, whereas the number of flavors with equal mass, thought of as the power of the 
fermion determinant is a continuous parameter. In the low energy partition function it is 
just the other way round. 



3 Supersymmetric Calculation of the Resolvent 



In this section we calculate the super-integrals in (|34D to obtain analytical expression 
for the resolvent and the microscopic spectral density. For integer u, this calculation 



was also presented in 1 11, 12]. To simplify the explicit calculations we will choose the 
parameterization 

f/^fC "lexpf" "V ,32) 



In this case the Haar measure is simply given by 

dU = d9dsdad(3. (33) 
This results in the partition function 

Z{z, z + J) = — ds dedadpe'"^-"' exp StrM 
27r J-oo Jcc 




(34) 



Where z and J are now microscopic variables, i.e. they are expressed in units of I/V4S0. 
The integration over s is over the complete real axis. For integer u the integration 
over 6 is over the interval [— 7r,7r]. For non- integer u the translational invariance of the 
^-integral is lost. It is recovered by extending the integration contour to include the 
intervals (— vr + ioo, — vr] and [n, vr + ioo). A picture of this integration contour, denoted 
by Cc, is shown in Fig. 1. 
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-71 + ioo 



71 + ico 



-71 



c. 



n 



Figure 1: Contour for the 6 integration. 



After performing the Grassmann integrations the partition function reduces to 

ZTT J-oo JCc 

= {z + J)K,{z)h^^{z + J) + zK,+^{z)h{z + J). (35) 

The normahzation of the partition function according to Z{z, J = = 1 follows from 
the Wronskian identity zKi,{z)I^^i{z) + zKi,^i{z)I„{z) = 1. The resolvent obtained by 
diferentiation with respect to J is, after using some identities for Bessel functions, given 
by 

E{z) = z{K{z)I,{z) + K_^{z)I,+,{z)) + ^. (36) 

The microscopic spectral density then follows from the discontinuity across the imaginary 
axis according to ([T0|) , 

Ps{x) = ^[E{ix + e) + T.{-ix + e)] 

ZTT 

= ^{Ju{xf-Ju+i{x)Ju-i{x)) + iy6{x). (37) 

The last term is the contribution from the u zero modes. For integer values of u the 
integration contour Cc can be replaced by the segment [— 7r,7r]. The restriction of the 
integration over 6 to this segment for non-integer values of u would have resulted in the 
wrong answer. This can be seen from the following representation of modified Bessel 
functions, 

1 /'^'^ sin UTT 
I^(z) = — e''°'^e'''^de-- / e-"™^^"-"^(is. (38) 

271 Jo 71 Jo 



11 



The microscopic spectral density and the expression for Ti{z) were first obtained from 



chiral Random Matrix Theory by means of the orthogonal polynomial method [31, 32 



and can also be obtained by means of the supersymmetric method I n that case 

the derivation, starting from the joint probability distribution of the eigenvalues, is also 
correct for non- integer values of u and is given by the expressions (|36[) and (p7|). 
The asymptotic expansion of the imaginary part of ^{z) is given by 



-2z 



2z 



+ 



-i{Aiy^ - l)(-l)^e-2^ 



z(-l)'^e-2-(4z/2- l)(4z/2-9) 



64^3 



+ 



i(4z/2 _ I)(4zy2 _ 9)(i9 _ 4z/2)(-l)^e 



-2z 



3 X 256^4 

and for the asymptotic expansion of the spectral density we find 



(39) 



1 

IT 



1 



cos(2a; 



TTZ/ 



2x 

2 



+ (4z/' - l)(4z/' -9) 



+ 
-6 



1 



4^(1 - sin(2x - TTZ/)) + (4z/2 - i)(4z/2 _ 9)-°'^^"^ 



TTU] 



64x3 



;i9 - 4z/2) sin(2x - ttu)) 
x4273! 



+ 



(40) 



The asymptotic expansion of the partition function, the resolvent, and the spectral density 
terminates for half-integer values of u. For example, for u = 1/2 only one oscillating term 
in (|40|) is nonvanishing suggesting that it can be obtained from a leading order saddle 
point approximation. 

Since for half-integer u the integral is localized on the critical points, we expect that 
the asymptotic expansion generated by the replica trick reproduces the exact answer. For 
other values of u we expect that the replica trick reproduces the asymptotic series to all 
orders in 1/x. 



4 Spectral Density via Fermionic Replicas 

In this section, we analyze the low-energy partition in the quenched case {Nf = 0). For 
n replica flavors in the sector of topological charge u it is given by (see ([T^)) pOj] 

Zl"\x) = f dUdet''Ue^^'^^+^~'\ (41) 

JUeU{n) 

In order to take the replica limit, n — ^ 0, the n-dependence of (^1]) has to be known 
explicitly. The closed form of Z^^'^x), in terms of a determinant of modified Bessel 
functions, 

Z(")(a;) =det(4+,_,(x)), z,j = l,---,n, (42) 

does not provide us with an explicit n-dependence. The explicit n-dependence can only 
be obtained for the large- mass and the small- mass expansion of Z^"-* (x) . It was obtained 
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m 



2^ , ^ for an expansion about the rephca symmetric saddle point using the method 



of Virasoro constraints. Our new result is the asymptotic expansion of the microscopic 
spectral density using replica symmetry breaking a la Kamenev and Mezard [T^ which 
will be discussed the the second half of the next subsection. The small mass expansion of 
(^) was also considered in |2^, but failed at the order for which logarithmic terms 
enter in the expansion. In the second subsection we derive the lowest order logarithmic 
term for the case of zero topological charge. 



4.1 Large Mass Expansion 

The asymptotic expansion of the microscopic spectral density is obtained from the large 
mass expansion of the finite volume partition function. To this end we expand the par- 
tition function (|41|) in powers of 1/x by means of a saddle point approximation. By 
diagonalizing U it can be easily seen that the saddle points are given by unitary matrices 
with eigenvalues ±1, i.e. by unitary matrices satisfying U"^ = 1. The solutions of this 
equation are highly degenerate. They can be organized in + 1 classes, U = Ip, where Ip 
is a diagonal matrix with p elements —1 and n — p elements +1 (with < p < n). The 
integrand does not depend on the submanifold U{n)/U{n — p) x U{p) of U{n) and the 
integration over this coset has to be performed exactly resulting in its volume Vn,p- In 
terms of the parameterization 

U = IpUoVU^\ (43) 

with Uq G U{n)/U{n — p) x U{p) it is clear that the integration can be restricted to 
U{n — p) X U{p). Summing over all saddle points the partition function is given by 

n „ 

Z(")(a;) = yVnJ dVJ{V)det''Ve-^^'^'^^^+^'"^\ (44) 

^ ' JV(iU{n-p)xU{p) 

where Ki,o = 1 and for p 7^ 

I ) Yt S^' ^ (27r)^("-^)F,^ (45) 

The integration over V should be thought of a saddle-point integral of a formal expansion 
of V about the identity to all orders. Below we will make this explicit for the different types 
of saddlepoints. The total number of saddle-points in the class p is ■ This factor is 

included as combinatorial factor in Vn^p- The volume of the coset U{n)/U{n — p) x U{p) 
is given by the ratio y(^(p^) y(^\^„_p)) ■ With the volume of f/(A;) given by (27r)^^ nj=ii' 
we obtain the volume factor p. 

From the exact expression of the partition function ( ^if ) in terms of modified Bessel 
function given in (|42|) it is clear that the asymptotic series of Z^{x) terminates for half 
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integer u. Let us investigate the asymptotic expansion about the saddle point Ip in more 
detail. Because the total number of degrees of freedom in U{n — p) x U (p) is equal to 
{n — p)^ + p^, the expansion of Zl"\x) is of the form 

~ e^"-'^')" (-] (1 + C f-j). (46) 

This result is valid for arbitrary u. The asymptotic series of the Bessel function /^_,_ i (x) 
terminates at l/x^"*"^/^. The result for the maximum power in 1/x occuring the expansion 
of the determinant (^2]) is particularly simple for v = \ and is given by 

- • (47) 
X J 

Let us consider this case in more detail. Since, as we will see below, only the saddle points 
for p = or p = 1 contribute in the replica limit, we only discuss these values of p. We 
observe that the maximum power and the minimum power in the asymptotic series are 
equal for p = saddle point. The asymptotic series thus has only one term and we find 
that 

X .IX n2/2 

Z\%{x)] ^ e-i-] . (48) 



p=0 \x 

For p = 1 and larger half-integer values of z/, more terms contribute to the asymptotic 
series, but it still terminates. For example, for p = 1 and = \ one finds from (^) and 
(^) that the difference between the maximum and minimum power in the expansion in 
l/\/x is n — 1. Therefore, the asymptotic series for Z^{x) cannot contain more than n 
terms. From numerical examples for small values of n, one indeed finds that in this case 
the coefficients of all n possible terms are nonvanishing. For the same reason as in the 



case of the GUE two-point function [14], we expect that the replica trick will give the 
exact result both for z/ = | as well as larger half-integer values of u. 

Clearly, the expression (Q) makes only sense for positive integer values of n . In 
order to analytically continue it we follow the work of Kamenev and Mezard [jl3l . They 
analytically continued the factorials in such that vanishes for p > n + 1. Then the 
sum in (|4^ ) can be extended up to infinity and the replica limit n — can be taken term 
by term in the sum over p. One finds 

hm FP ~ nP, (49) 



SO that only the terms p = and p = 1 of the, to infinity continued, sum in ( ^^ survive. 
In these two cases we obtain F^^q = 1 , F^^q = n. Notice, that the continuation of the 
sum over p to infinity explicitly breaks the replica symmetry p n — p. Of course, the 
group integral in (30) must be also continued to non-integer values of n. 
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The previous discussion suggests the definition 

Z^-\x) ^ [Z^'\x))^^^ + n [Z^J'Kx))^^^ . (50) 

We will first consider the contribution for p = which originates from expanding U 
around the identity matrix Jq. Although not necessary, it turned out to be convenient to 
parameterize Vp=o ^ U (n) according to 

1 + iH/2 , , 

where H is an Hermitian n ^ n matrix. From a diagonal representation of lp=o one can 



easily show that 



In the replica limit, the Jacobian can be ignored and one simply has dVp=o = dH. The 
p = contribution to the partition function is thus given by 



(53) 



In order to keep track of the powers of x it is convenient to rescale H — > H/ a/x which leads 
to another Jacobian x~'^^^'^ that also vanishes in the replica limit. In terms of x = mEoV^, 
the mass dependence of the condensate is given by 



Sp=o(x) ,.119 



lim 



So n-on^(")(a;)9x 



so that to order 1/x^ we need to collect terms to order 1/x^ in the expansion of the 
partition function. Using the expressions 

(det\/p=o)" = e'^'^e^V i^+' V . (55) 
we find the result (only terms up to order 1/x^ are displayed) 



V " /p=o ^4 J 



X 2 

:i - ^(Trff)^ + ^(Trff)(Tri7)^ + ^(M)^) 

'1 + — Trif^ ^Tri/^) 

8x 32x2 > 

-) e"(l-„- + (2„3 + „)-) + 0(^). (56) 
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The Gaussian integrals have been calculated using the trace correlators given in the Ap- 
pendix. To this order, the Jacobian in ( ^21) gives rise to an extra factor (1 — n^x/A). With 
inclusion of this terms the 1/x corrections vanish for p = \. All terms of order l/x^ are 
at least of order n^. The resolvent obtained from (|5^ 



Sp=o(a;) 



1 + 



(4z/2-l)(4z/2-9) 



+ 



(57) 



Eo 8x2 128x4 

agrees with the large mass expansion of the quenched condensate obtained via other 
methods [^. We observe a cancellation of the odd powers in 1/x. This is in agreement 
with the analytical expression for E(x) given in (^) which can be rewritten as 



Sfx 



X 4z^ 

-- = -((2 - —)I^{x)K,{x) + /,+i(x)ir,+i(x) + /,_i(x)ir,_i(x)). 
Lo 2 x^ 



(58) 



From the asymptotic behavior of the Bessel functions, one easily derives that the asymp- 
totic series of this expression is an expansion in powers of l/x^. Also notice the cancella- 
tion of the v/x term in the large— x asymptotic expansion. 

Next we consider the more subtle contribution of the saddle point given by the diagonal 
matrix /i with one element equal to —1 and n~l elements equal to 1. A parameterization 
of U{n — 1) X f/(l) that is convenient for the expansion about the saddle point is given 
by 



l+iH/2 
l-iH/2 



l+ih/2 
l-ih/2 



(59) 



where, now, H is a hermitian (n — 1) x (n — 1) matrix and h is a real variable. Because 
Uq e U{n)/U{n - 1) X f/(l) we have that 

U = hUoVp=iU^' = UohVp=iU^\ (60) 

The measure, in terms of the coordinates H and h, can be obtained by diagonalizing 
U and H with unitary transformations Ui and W, respectively. If the eigenvalues are 
denoted by e*^ and hk, in this order, we find the measure (notice the plus sign in the last 
factor) 



dU 



n 

k<l 



Ilk<i \ hk — hi\ 



'l + hy4)Uk{l + hl/A) 



The last factor can be written as 



n—l 



n 



l + ihk/2 l+ih/2 



ihk/2 



ih/2 



dh Yl dhk dUodV. 

k 



n 

k 



l + ihk/2 ^ l + ih/2 



ihk/2 l-ih/2 



i n— 1 



det^(l + hH) 



(l + /iV4)"-idet(l + iJV4) 



(61) 



(62) 
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Using that Y[k<i \hk — hi\'^ Yik dhk dV — dH we thus find the measure 



dU 



tn—l 



dUo- 



dhdH det\l + hH) 
l + /iV4)Met"(l +ifV4) 



= J{H, h)dhdHdUo. 



(63) 



The integrand does not depend on Uq and the integration over these variables just gives 
the vohimc of the coset which, together with the combinatorial factor, combines into 
the factor \/"„ i discussed in the first part of this subsection. The p = 1 contribution to 
the partition function is thus given by (with one overall minus sign from orienting the 
/i-integration from — oo to oo) 



X 



-{-IY{StiY-^ j dHdh — 



det2(l + hH) 



X 



(l + /iV4)Met"(l + i/74) 



1 - ih/2 



1 - iH/2 



(64) 



Since the factors in the denominator of the Jacobian can be ignored in the replica limit 
only the following terms in the expansion of the Jacobian contribute to order 



J = 4 



n-l 



1 + 



2 ^16 



2(Tr//)2 - Tr//^ 



+ 



(65) 



It is instructive to perform the calculation to leading order in 1/x. In this case the 
following terms should be collected 



{Z^''\x))^^^ = -(-l)"(87r)"-^e("-')^+'^'^y'rfi/rf/ie-tT''^' 



X 



1 + flV//^ - ^/i^ - ^((Tr//)2 + /i^ + 2ihT,H) + 
8 8 2 2 



(66) 



The last two terms vanish upon integration. Thus, the Jacobian contributes only at the 
next to the leading order. The saddle-point integrations can be performed conveniently 
by rescahng h and H according to /i — > h/\/^ ; H — > H/^Jx. This results in 



i(-l)^(87r)"- ( 



n-l / ^" \ {n-2)x+vn 



X 



1 



1 + (2(n - 1)^ + (n - 1))- - - - n{n - 2) 



'2x 



(67) 



Now we are in a position to calculate the contribution of the p = 1 saddle point to the 
chiral condensate. 



Sp=i(x) 1 1 d . 

= lim— — — — — n 

Eo n^o n Zl}{x) dx^ 



n). 



X 



(68) 
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Using trace correlators given in the Appendix, the expansion in 1/x can be easily extended 
to order Xjx^. For the mass dependence of the chiral condensate we obtain 

—-^ = hm (x)—Zl^(x) 

^ i(-l)"e-^" ^ (4t/'-l)(l-z(-l)"e-'") z(-l)^e-2-(4z/2-l)(4z/2-9) 



2x 64x^ 
(4z/2 - I)(4zy2 _ 9) [z(-l)-e-2-(19 - 4i/2) _ 6] 

3 X 256x4 ■ 



(69) 



All terms ~ ie~^^ originate from p = 1 saddle point. Before calculating the spectral 
density by taking the discontinuity of S(a;) we wish to point out that the expression (|69|) 
has been obtained under the assumption that Re(a;) > 0, so that we cannot calculate 
the discontinuity from the difference of S(iA + e) and S(iA — e). The reason is that for 
x = 'iA + ethe dominant saddle-point is given by Jq and for x = iA — e it is given by 
—Iq. The infinitesimal increment thus breaks the replica symmetry between Iq and —Iq. 
The sum over p in P^ ) has been extended to oo consistent with the breaking of the 



replica symmetry p ^ n — p hj the saddle point Iq. Thus, the infinitesimal increment 
is necessary to resolve the ambiguity as was also the case in the original calculation of 



T3[| . The expression for S(iA — e) can be obtained from a replica-symmetry breaking 



solution for which —Iq dominates. The final result for for the asymptotic expansion 



of Ps(A) coincides with (^) obtained from the expansion of the analytical result 
(A/2)(J2(A) - J,+i(A)J,„i(A)) up to l/A^ for A > 0. 

We observe that it requires a great deal of effort to derive the asymptotic expansion of 
the oscillating contribution to the spectral density by means of the replica trick. This is 
especially true due to the lack of the Virasoro constraints for z/ 7^ 0. Those constraints are 
a key tool p4[ in simplifying the calculations for the mass expansions of the condensate in 
the sector of vanishing topological charge {u = 0). The simplifying feature in the sector 
with vanishing topological charge that the partition function can be shown to belong 
to the universality class (see [Q) of the generalized Kontsevich model with potential 
V(x) = and satisfies the same Virsoro constraints. For z/ 7^ the one-link integral 
depends also on det J and det which is also the case when we have integrals over 
SU{Nf) instead of U{Nf). Therefore, a possible generalization of the Virasoro constraints 
to include the case u ^ would be certainly welcome not only to reduce our efforts in 
perturbative calculations but also as an identification of the universality class of the QCD 
finite volume partition function at nonzero topological charge. 

4.2 Small Mass Expansion 

The situation for the small mass is much more complicated. Before discussing the com- 
plications we recall that the partition function Z^{x) has been extensively studied in the 
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context of lattice QCD where it is known as the one-link integral (see for example |^ for 
an updated review of the subject). In that context one considers unitary matrix integrals 
of the form 

Z(J,jt)=/ rf[/e^K^^'+^^^), (70) 

JueU{n) 

where J is a general nxn matrix. The partition function with such potential is a function 
of the eigenvalues of the matrix JJ^. The one-link model exhibits two phases according 



to p 



— j= < 1 (weak coupling), 
fc=i 2v Afc 

n ^ 

— y= > 1 (strong coupling). 



(71) 



In our case, with partition function given by (^21), the eigenvalues of J^J are given by 
From the above we see that, for x oo, we expect to take the replica limit, 
n 0, and remain in the weak coupling regime while, for s — > 0, it is not clear whether 
the replica limit can be taken without crossing a phase boundary. These problems are 
reflected in logarithmic singularities of the small x expansion of the valence quark mass 
dependence of the partition function and its derivatives. In the replica limit of the 
expansion coefficients could be derived up to the order for which terms of the form x^ log x 
are absent. These singular terms could be related to the presence of de Wit-'t Hooft 



poles 



Below we consider the replica limit of Z'^{x) for v = Q. In that case logarithmic terms 
already enter to lowest order in the expansion. We thus consider the small x expansion 
of the partition function Z"^q(x) = det(/j_j(x)). By expanding the Bessel functions we 
obtain 

n / 2 \ 1 oo 
fc=0 \ ^ / k=n+l 

The first sum can be recognized as an incomplete exponential, but we were not able to 
determine the coefficients Ck^n in general. In order to expose the n-dependence of the first 
term, we rewrite the incomplete exponential as an incomplete gamma function 




fc=0 



e^r(n + l,a;V4) 

r(n + i) ■ 



(73) 



From the small mass expansion of the incomplete gamma function 



' 9\ k+2n 
X ^ 



r(„ + 1. .74) - dt re- ^ r(„ + 1) + g j (74) 
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and neglecting terms which will vanish in the replica limit, i.e., 



we find 



4 



1 + n 



r(n + l) 

2\ " 1 



n + k 



1 - 7n + 0{n^) 
1 /2, /i 1 



0(n2), 



(75) 



:i-e^'^/^)(21og| + 7) + e-^/^E — 

fc=i \ 



2 \ ^ 



k\k 



+ 



k=n+l 



2k 



(76) 



Because the coefficients Cfc^„ are unknown only the lowest order terms of the small mass 
expansion can be calculated, 

(77) 



X 



1 - riy logx + O(x^) + 0{n'^) 



For the replica limit of the condensate given by. 

1 



lim 



dZ^{x) 



onZnix) dx 



we thus obtain 



E(x) = —X \ogx + 0{x^). 



(78) 



(79) 



This is indeed the correct leading order term of the small mass expansion of x{Io{x)Ko{x) + 
Ii{x)Ki{x)). The linear behavior of the microscopic spectral density at the origin is 
reproduced by taking the discontinuity across the imaginary axis. 



PsiX) = ^ + OiX'). 



(80) 



We did not succeed to generalize this calculation to arbitrary u, but we expect that 
the logarithmic terms can be obtained in a similar fashion. In particular, the first n + u 
terms of the expansion seem to follow a simpler pattern than the coefficients of the higher 
powers. 



5 Bosonic Replicas 



In view of the seemingly different role [0, 14 1 played by bosonic and fermionic replicas and 
the fact that the supersymmetric method uses both compact and noncompact variables 
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it is natural to try to reproduce the results of last section by introducing additional n 
replicas of bosonic quarks of mass m instead of fermionic ones. In this case, the condensate 
for Nf = is given by 

= lim^-^^lnZ(-")(x) , (81) 
n-*o (-n) V^dm 

where Zl~"\x) is given in (^). We will show that, for large masses, bosonic replicas can 
be used to reproduce the asymptotic expansion of the chiral condensate but they fail 
to reproduce the microscopic spectral density for a subtle reason which we will explain 
below. 

It is convenient to express the bosonic partition function Z^~^\x) in terms of the 
eigenvalues of the Gl{n)/U{n) matrices. Up to an overall constant we have, 

= r fl dsklKe'" - e'^){e-'' - e-'Oe"'^^'^=i '°'^''+''^^=i (82) 

k=l k<j 

Notice in particular that the measure (including the Vandermonde determinant) is invari- 
ant under the symmetry Sk Sk + t which is a remnant of the Gl{n)/U{n) invariance. 
The fermionic partition function Z^\x) is given by a circular ensemble (for integer z/) 
0, 



= r n dOk n(e'^' - e^^O(e"'^' - e''^ ^)e'^"k=^''°'^ ^"^^^^K (83) 

-'-'^ k=l k<j 

In the noncompact case, the solutions of the saddle point equation sinhs^ = are 
given by Sk = 0, iivr, ±i27r ■ ■ -. Thus, in principle we might have a variety of saddle 
points which should be all taken into account in large x expansion. However, we will 
argue that only Sk = solution contributes to the large x behavior of Zj~^\x). Our 
discussion is based on the n = 1 integral where a steepest descent analysis can be easily 
carried out. In this case, the bosonic partition function is given by 

roc 

Zl-^\x)= dse-'"'°'^^'^+''' = 2K,{x) (84) 



From the asymptotic expansion of Ky{x) it is clear that only the saddle-point at s = 
contributes to the integral. The spectral density can be calculated from the resolvent at 
X = ±iX + e so that the integral above can be identified with the modified Bessel function 
KJx). 



On the other hand, for one fermionic replica the partition function is given by (assum- 
ing integer u) 

Zi'\x) = 277/, (x). (85) 
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Im s 




Figure 2: Steepest descent curve passing through the saddle point s = for A < (sohd) 
and for A > (dashed), where we assumed x = i\ + e. 



In this case two saddle-points contribute to the asymptotic expansion, one ~ e^, and the 
other one ~ e~^, which are both essential to recover the oscillatory contributions to the 
microscopic spectral density. 

Let us analyze in detail the saddle-point calculation of the bosonic partition function 
for one replica. The stationary phase condition that the imaginary part of the action, 
— iAcoshs, is constant results in the following curve in the complex s-plane through the 
saddle point at s = 0. Clearly, the integration contour (the real axis) can be deformed 
into the steepest descent curve of fig. 2 (depending on the sign of A) and no other saddle 
points need to be considered. A similar analysis for the compact case shows that both 
the saddle points at ^ = and at = vr have to be taken into account. 

A similar situation arises in the correction to the semicircle law of the Gaussian 
Unitary Ensemble (GUE) of hermitian matrices. Inside the semicircle the saddle point 
solutions of the n = 1 fermionic replica are horizontally aligned in the complex plane, 
while (see also a comment in [|I4|) the solutions for the one bosonic replica are vertically 
aligned and only the p = saddle point contributes. If we insist on taking into account 
all vertically aligned saddle points we get an incorrect result. However, we have checked 
that outside the semicircle the situation is reversed (see figures in |^ for the fermionic 
case). In this case the bosonic replicas correctly reproduce the semicircle law and its 
leading exponentially decreasing correction outside the semicircle (see also [^) while 
the fermionic replica calculation gives an incorrect result. In the case of the microscopic 
spectral density analyzed in this work, we are always inside the semicircle even in the 
limit X ^ oo which explains, assuming the same pattern of as for the GUE, the failure of 



22 



the bosonic replica calculation in reproducing the oscillating part of the spectral density. 

If we are just interested in the mass dependence of the chiral condensate it is sufficient 
to only take into account the p = saddle point and, one can convince oneself that 
bosonic and fermionic replicas produce the same large x result as follows. The integrand of 
the non-compact partition function is obtained by transforming the integration variables 
according to 6^ = isk and replacing x —>■ —x. For large x, the saddle-point of both 
partition functions is at 6*^ = and Sk = 0, respectively, and one finds that —X) 
and Zl~^\x) have the same asymptotic expansion for the chiral condensate (see (p7|)). 
Clearly this proof is formal because it assumes that the replica limit can be interchanged 
with the operation x —x. However, we have explicitly checked to the order 1/x^ that 
the noncompact partition function Zl~"'\x) leads to the same asymptotic expansion. 

6 Conclusions 

We have investigated the replica trick for the microscopic spectral density of the QCD 
Dirac operator in the quenched limit. The advantage of working with fermionic replicas 
is that this theory corresponds to QCD with n fiavors of equal mass m. Because the 
low energy properties of this theory have are well understood, the starting point of this 
approach has a firm basis. The valence quark mass dependence of the chiral condensate 
and the spectral density of the QCD Dirac operator, however, are only obtained in the 
limit n — i> 0. The existence of this limit has been debated for many years and the 
investigation of its nature has been the main topic of this article. 

The alternative approach to obtain the QCD Dirac spectrum is the super symmetric 
method. Although in principle rigorous, one might raise the question whether this theory 
with bosonic ghost quarks might have unusual properties as for example the spontaneous 
breaking of supersymmetry. Our results show that this is not the case. The low-energy 
limit of this partition function is completely dictated by chiral supersymmetry. The power 
of the supersymmetric method is that one obtains rigorous nonperturbative results such 
as, for example, the spectral density in the microscopic region. 

The replica trick, on the other hand, requires an explicit n dependence, and, up to now, 
only perturbative results have been obtained. Exact results have only been derived in cases 
where the perturbative series consists of only a finite number of terms. Our results for 
the microscopic spectral density confirm that the asymptotic series of its non-oscillatory 
part can be obtained from an expansion about the replica symmetric saddle-point. Both 
fermionic and bosonic replicas give the same result and are in complete agreement with the 
asymptotic expansion of the exact result. The same is true for the asymptotic expansion 
of the resolvent away from the imaginary axis (where the eigenvalues are located). Things 
are different for the the oscillatory part of the microscopic spectral density. In this case 
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the correct asymptotic expansion is obtained only if a saddle point that breaks the replica 
symmetry is taken into account. This additional saddle point only exists for fermionic 
replicas. For bosonic replicas, only the replica symmetric saddle point contributes in the 
saddle point calculation, and therefore this approach does not reproduce the asymptotic 
expansion of the oscillatory part of the spectral density. A similar observation has been 
made for the application of the replica trick to the Wigner-Dyson ensembles. 

In the asymptotic expansion of the supersymmetric partition function also two saddle- 
points have to be taken into account. In the boson-boson component of the Goldstone 
manifold only one saddle point contributes, but as is the case for one fermionic replica, 
we have to take into account two saddle-points in the fermion-fermion component of the 
Goldstone manifold. 

Is it possible to go beyond perturbation theory using the replica trick? One indication 
in favor of an affirmative answer to this question is that we have reproduced the leading 
order logarithmic singularity of the small mass expansion of the resolvent. Except for 
these logarithmic terms, the small mass expansion of the resolvent is a convergent series 
which may be summed to obtain the exact result. The large mass expansion, on the other 
hand, is an asymptotic series which cannot be summed and cannot provide us with the 
exact result. Of course, for half integer v when the asymptotic series terminates, an exact 
result is obtained from the replica trick. 

The exact answer for the resolvent in the microscopic region shows a compact- noncompact 
dichotomy. This dichotomy is natural in the supersymmetric approach where the com- 
pact part of the Goldstone manifold is associated with the fermion-fermion sector and 
the non-compact part is associated with the boson-boson sector. In the fermionic replica 
trick this dichotomy is not at all clear but might be hidden in the n — > limit which is 
given by an integration over a 1 x 1 matrix and and an integral over an {n — 1) x (n — 1) 
matrix. This might be another hint that it is possible to go beyond perturbation theory 
within the replica framework. 

Finally, we hope to have convinced the reader that the supersymmetric method is 
the only (j— model approach that can provide us with rigorous exact results. Even the 
calculation of a small number of terms in the asymptotic expansion within the replica 
approach requires a tremendous effort in the case of broken replica symmetry. 
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Appendix 

We are interested in expectation values of traces of powers of nxn Hermitian matrices 
with matrix elements distributed according to the Gaussian Unitary Ensemble (GUE). 
Such averages are given by 

J an e 2 

IZom=i (rfa;.e-H)n,<,(x.-x,)2 
= (86) 

where ui = TtHK All correlators Qg{pi, ■ ■ ■ ,pg) can be calculated recursively starting 
from fii(O) = n. The recursion relations are derived from integrals of total derivatives 
(Schwinger-Dyson equations), 

n n I n n \ 

n rfx, E dk xl^' E ■ ■ ■ E K = 0. (87) 

-°°i=l k=l \ h=l jt=l ) 

It is useful to keep in mind the identity 

E5.(^r'Ay =A^^X:^^^-^ (88) 



fc=i i=\ 



As a sample calculation let us derive ^23(1, 1,4). Choosing 6 = 2,ri = 4, r2 = 1 and 
a = — 1 in (B7D we find 



r]3(l, 1,4) = 4^2(1,3) +nfii(4). (89) 
Choosing now 6 = 1, a = — 1 and ri = 3 or ri = 1, respectively, we get 

^2(1, 3) = 3fii(2), 

^2(1,1) = n. (90) 
Finally, from 6 = and a = 2 or a = 0, respectively, we deduce 

fii(2) = n^, 

fii(4) = 2r2r]i(2) + ^2(1,1). (91) 
This system of equations results in 

^3(1, 1,4) = 2n^ + 13n2. (92) 
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For some special correlators we can easily derive a general formula, e.g., choosing b = 
2k, a = —1 and ri = r2 = ■ ■ ■ = r2fc = 1 we find 



n2fc(l,l,---,l) = n''{2k-l)\\. (93) 

Clearly all figipi, ■ ■ ■ ,Pg) with J2iPi being odd vanish identically. Besides (0) we have 
used the following moments in the calculation of this paper. 



ill 
ill 

ill 
Oi 

^2 

VLi 

^3 

^3 
^2 
^4 
^^3 



:0) = n, 
2) = n\ 

4) = 2n^ + n, ^^2(2, 2) = n^ + 2n'^, 

1,3) = 3^2, 1]3(1, 1,2) = n^ + 2n, 

-Q) = 5n^ + 10^2, 1^2(2, 4) = 2n^ + 9^^ + 4n, 

3.3) = 12^3 + 3n, ^2(1,5) = lOn^ + 5n, 
1,1,4) = 2n^ + 13n2, ^4(1, 1, 1, 3) = 9n^ + 6n, 

8) = Un^ + 70n^ + 21n, 1^2(2, 6) = 5n^ + 40n^ + 60n^, 

4. 4) = An^ + 40n^ + Qln^, Ogjl, 1, 6) = Sn^ + 70^^ + 30n, 

1, 3, 4) = Qn^ + 75n^ + 24n, (^5(1, 1, 1, 1, 4) = Qn^ + 75n^ + 24n, 

2, 4, 4) = + 72^6 + 381n^ + 488n2, 

4, 6) = lOn^ + 169n5 + QlOn^ + 156ri, 

1, 1, 4, 4) = 4n^ + 88n^ + 661n^ + 192n, 

A, 4, 4) = 8n^ + 228n^ + 220271^ + 6517n^ + 1440n. 



The sum of the coefficients of the correlators can always be checked by means of the n = 1 
case where all correlators reduce to one dimensional Gaussian integrals. 



(Ep. - 1): 



The coefficient of the highest order power in n of the correlators ^1(2^) can be checked 
as follows. First, in the limit n —>■ 00 one can solve the loop equation (or Virasoro 
constraints) which, for our Gaussian potential, gives the equation 

~[p-Xi 2 

where p is a positive number assumed to be large (outside the semi-circle). After expand- 
ing both sides of this equation in powers of 1/p we find 

lim VLii2k) = ^ — '—^ — 97 
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